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Abstract 

The idea of working with a near-critical phase-separated liquid mixture whereby the surface 
tension becomes weak, has recently made the field of laser manipulation of liquid interfaces a 
much more convenient tool in practice. The deformation of interfaces may become as large as 
several tenths of micrometers, even with the use of conventional laser power. This circumstance 
necessitates the use of nonlinear geometrical theory for the description of surface deformations. 
The present paper works out such a theory, for the surface deformation under conditions of axial 
symmetry and stationarity. Good agreement is found with the experimental results of Casner and 
Delville [A. Casner and J. P. Delville, Phys. Rev. Lett. 87, 054503 (2001); Opt. Lett. 26, 1418 
(2001); Phys. Rev. Lett. 90, 144503 (2003)], in the case of moderate power or a broad laser beam. 
In the case of large power and a narrow beam, corresponding to surface deformations of about 50 
micrometers or higher, the theory is found to over-predict the deformation. Possible explanations 
of this discrepancy are discussed. 

PACS numbers: 42.25.Gy, 42.50.Vk, 82.70.Kj 
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I. INTRODUCTION 



The formation, deformation, and breakup of fluid interfaces are ubiquitous phenomena 
in nature Q . One special group of effects which implies the so-called finite time singularity 
Q], has as its most common example the breakup of pendant drops driven by gravity ^, Q|. 
If external fields such as electric or magnetic fields are present, one has in addition extra 
control parameters. Thus we may mention that interface instabilities driven by electric fields 



are important for many processes such as electro-spraying [7|, ink-jet printing [8|, or 
surface-relief patterning [9j]. A uniform magnetic field can also be useful, for instance for 
the purpose of of forming elongated magnetic droplets |lol |. These deformations, as well as 



11 



have been used 



13, 



those induced by the acoustic radiation pressure on liquid surfaces 
to explore the mechanical properties of fluid interfaces in a non-contact way 

It is noteworthy that laser-induced deformations of the interfaces of soft materials have 
not received the same amount of interest in the past. Most attention has been given to test- 
particle global effects such as optical levitation and trapping - cf., for instance, Refs. [3, 



16(. The reason for this circumstance is simple: deformations of fluid interfaces by optical 



radiation are ordinarily weak. For instance, in the classic experiment of Ashkin and Dziedzic 



17| . a pulsed single transverse mode doubled Nd:YAG laser (20 pulses per second) was 
focused vertically from above on a water surface. The wavelength of the incident wave was 
A = 0.53 fim, the peak power was P ma x = 3 kW (low enough to make nonlinear effects 
negligible), and the duration of each pulse was 60 ns. The beam radius at the waist was 
reported to be very small, u Q = 2.1 /im, but has most likely have been somewhat larger, 
Uq = 4.5 /im; cf. the dicussion in Ref. |l8j . The rise of the water surface was small, about 
0.9 (Urn, occurring at t ~ 450 ns after the onset of the pulse. The physical reason for this small 
surface elevation is evidently the large surface tension a =0.073 N/m between air and water. 



[The theory of the Ashkin-Dziedzic experiment was worked out in Refs. 18|, ll9j].] If we go 
to the more recent experiment of Sakai et al. j^], we will find that the surface displacement 
was even more minute. In this case the displacement was induced by a continuous wave 
(cw) Ar + pump laser (wavelength in vacuum Ao = 0.514 /xm, maximum power P max =0.5 
W), and was probed with a He-Ne laser. For a beam power P= 0.3 W and a beam waist 
ujq = 142 fim the elevation of the surface was found to be extremely small, about 2 nm. 
Generally speaking it is of considerable interest to be able to probe the displacement of 
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fluid interfaces in a way that is non-contact, i.e., which avoids a direct touch of the fluid 
by mechanical devices. The only acting force on the fluid surface is thus the radiation 



force. As discussed in Ref. 



2lj . for instance, this kind of force can measure locally the 



221 or 



micromechanical properties of soft biological systems because artificial membranes 
cells |23( can be highly deformable. Sizable effects of the radiation pressure should thus 
be measured, facilitating the characterization of surface elasticity properties. A major step 
forward was to recognize that one can reduce the surface tension considerably by working 
with a phase-separated liquid mixture, close to the critical point. In this way "giant" 



e have shown t 
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hat 



24, 



25 
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deformations can be achieved. Recent experiments of Casner and Delvil 
the displacements can in this way reach several tenths of micrometers 
(cf. also the recent review [28]). Also, there are seen in the experiments rather remarkable 
asymmetries with respect to the upward/downward direction of the laser beam 127 1. 

The giant deformations make it necessary to use nonlinear theory for the description of 
the surface deflection. The purpose of the present paper is to present such a calculation, for 
the typical case where the elevation is axially symmetric (a condition almost always satisfied 
in practice), and in addition stationary. Comparison with the mentioned experiments will 
be made. We shall moreover assume that the laser beam is incident from below only. The 
mentioned up/down asymmetry will thus not be treated. 

For completeness we emphasize that we are considering nonlinear theory only in a geo- 
metrical sense, in order to describe the interface bulge. The electromagnetic theory as such 
is kept on the conventional linear level, as we are only considering moderate laser intensities. 



The theory for linear deformations has been worked out before 
solution for the deflection of the surface is reproduced in Eq. (JT7j) . 



30, 2i, 3- 



The analytic 



II. DERIVATION OF THE GOVERNING EQUATION 



A. Basic electromagnetic theory 



We begin by writing down the expression for the electromagnetic volume force density f 
in an isotropic, nonconducting and nonmagnetic medium [if! [ffi]: 

k — Id 



c 2 dt 



E x H . 



(1) 
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Here E and H are the electric and magnetic fields, p is the mass density of the medium 
(fluid), e is the permittivity, and k = s/sq is the relative permittivity where Sq denotes the 
permittivity of vacuum. 

Let us comment on the various terms in Eq. (jT|). beginning with the last term. This 
term is called the Abraham term, since it follows from Abraham's electromagnetic energy- 
momentum tensor. The term is experimentally detectable under special circumstances at 
low frequencies [l^, but not at optical frequencies, at least not under usual stationary 
conditions. The Abraham term simply fluctuates out. 

The middle term in Eq. is the electrostriction term. When seen from within the 
optically denser medium (the medium with the highest n), the electrostriction force is always 
compressive. Whether this kind of force is detectable in a static or a stationary case, 
depends on whether the experiment is able to measure local pressure distributions within 
the compressed region or not. Moreover, in a dynamic case the velocity of sound is an 
important factor. If the elastic pressure in the fluid has sufficient time to build up, then 
the electrostriction force will not be detectable when measuring the gross behavior of a fluid 
such as the elevation of its surface. Such is usually the case in optics. The time required 
for the counterbalance to take place, is of the same order of magnitude as the time needed 
for sound waves to traverse the cross section of the laser beam. For a beam width around 
10 /mi, this yields a time scale for counterbalance of the order of 10 ns. For instance, in the 
Ashkin-Dziedzic experiment Il7f a detailed calculation verifies this time scale; cf. Fig. 9 in 

Q. 

Another point worth mentioning in connection with the electrostriction term is that that 
we have written (de/dp)s as an adiabatic partial derivative. This seems most natural in 
optical problems in view of the rapid variations of the field, at least in connection with laser 
pulses. In many cases it is however legitimate to assume that the medium is nonpolar, so 
that we need not distinguish between adiabatic and isothermal derivatives. The permittivity 
depends on the mass density only. Then derivative can be written simply as de/dp, and is 
calculable from the Clausius-Mossotti relation. In this way we can write Eq. in the 
following form, when omitting the last term, 

f = ~s E 2 Vk + if-oV [E\k - 1)(« + 2)] . (2) 

Finally, we have the first term in Eq. which may be called the Abraham-Minkowski 



force, since it follows equally well from the Abraham and the Minkowski energy-momentum 
tensors: 

f AM = -\e*E*VK. (3) 

This is the only term that we have to take into account in practice in optics, under usual 
circumstances. We see that this force is equal to zero in the homogeneous interior of the 
medium, and acts in the inhomogeneous boundary region only. By integrating the normal 
component of the Abraham- Minkowski force density across the boundary, we obtain the sur- 
face force density which can alternatively be evaluated as the jump of the normal component 
of the electromagnetic Maxwell stress tensor. 

In the following we focus the attention on the force term in Eq. (JHJ) only. 



B. Surface tension and radiation forces on a curved surface 

Let us assume that there is established a stationary curved surface z = h(x, y) distin- 
guishing two fluids, a lower fluid (1) and an upper fluid (2), the equilibrium position being 
determined by the balance of gravity, surface tension, and radiation pressure. The undis- 
turbed position of the surface is the xy plane. Because of the surface tension coefficient a, 
there will be a normal stress proportional to the mean curvature of the surface: 

P2 _ Pl = <7 (_L + _L), (4) 

Ri and R2 being the principal radii of curvature at the surface point considered. If R\ and 
i?2 are positive, P2~ Pi> 0. This means that the pressure is greater in the medium whose 
surface is convex. It is useful to have in mind the following general formula for the mean 
curvature 1/R\ + l/R^- 

1 1 h xx (l + h^) + h yy (l + h 2 x ) - 2h xy h x h y 

R' 1 + R'2~ (1 + ^ + ^2)3/2 > U 

where h x = dh/dx, etc. Our convention is such that the curvature is positive if the surface 
is concave upwards. Also, we note that the unit normal vector n to the surface is 

n=(l + hl + hl)- 1 ' 2 (-h 9 ,-h v ,l). (6) 
The normal points upwards, from medium 1 to medium 2. 
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Assume now that there is a monochromatic electromagnetic wave with electric field vector 
EW(r)e -l£J * incident from below, in the positive z direction. The direction of the incident 
wave vector kj is thus given by the unit vector 

ki = (0,0,1) (7) 

in medium 1. When this wave impinges upon the surface, it becomes separated into a 
transmitted wave E^ and a reflected wave E^, propagating in the directions of k t and k r , 
respectively. We assume, in conformity with usual practice, that the waves can locally be 
regarded as plane waves and that the surface can locally be regarded as plane. The plane of 
incidence is formed by the vectors kj and n; we call the angle of incidence 9i and the angle 
of transmission Q t . See Fig. 1. Moreover, we let En and Ej_ be the components of E parallel 
and perpendicular to the plane of incidence, respectively. The expressions for the energy 
flux transmission coefficients Tji and Tj_ for a plane wave incident upon a boundary surface 



are (cf. 



p. 496): 



n 2 cos#t / -^n \ sin 26 i sin 26 t 



2 



11 m cosft ^ Bf ) ~ sin 2 (^ + 6 t ) cos 2 (^ - 6 t ) ' 1 ' 

= n2 cos^t / E^_\ 2 = sin2^sin2fl f 
± ni cos0« ysf J = sm 2 (6 t + 6 t ) ' [ ' 

When dealing with an unpolarized radiation field, one usually averages over the two polar- 
izations and represents the transmission coefficient by the single entity 

(T) = i(Tn+T ± ). (10) 

Consider now the electromagnetic surface force density, which we will call II. As mentioned 
above, II can be found by integrating the normal component of the volume force density 
across the surface boundary layer. From Eq. Q it follows that the surface force acts normal 
to the surface, and that it is directed towards the optically thinner medium. 
We introduce the intensity / of the incident beam, 

/ = en lC (E^ 2 } (11) 

(in the case of azimuthal symmetry / = I(r)), and let a denote the angle between E^ and 
the plane of incidence, 

Ef = E {i) cos a, Ef = E {i) sin a. (12) 



Then, we can write the surface force density as 

II = .l ^-rcj cosft r 2 + cog2 cog2 a + ^ g . n2 , n 
2c n 2 cos 6^ L 

When EW = E[| or EW = (i.e., a = or 7r/2) it is often convenient to express II as 



n 



— COS 2 0*1 1-4 /; 



tan#,- 



T n, 



(14) 



where R 



<■ \ tan 6^ 

1 — T is the reflection coefficient. This expression has been derived before 

lydrodynamic nonlinear case. In connection with the 



H 27, 0- 

It holds also in the 
mentioned Bordeaux experiments 
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24, 



25, 



2d El 

, the upper liquid was always the 



optically denser one. Thus n% > n\, the direction of II is antiparallel to n, and the force 
acts downward, normal to the surface. 
The case of normal incidence yields 

4nin 2 



(n 2 + nif 



(15) 



2n 1 In 2 -n 1 
II = n. 

c n 2 + ni 



(16) 



C. Cylindrical symmetry 

We henceforth assume cylindrical symmetry, using standard cylinder coordinates (r, 9, z). 
There is no variation in the azimuthal direction, dh/d9 = 0. With the notation h r = dh/dr 
we have 



cos 9i 



1 



sin 9 s 



h r 



Together with analogous expressions for 9 t this can be inserted into Eq. (|T3j) to yield 

2nJ(r) 1 - a 



n 



where a is the relative refractive index, 



c 1 + a 



f{hr,a)n, 



and f(h r , a) is the function 



f(h r ,a) 



a = n 1 /n 2 < 1, 



[1 + a) 



(17) 



(18) 



(19) 
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x < sin a + 



1 + (3-a 2 )/i 2 + {2-a 2 )h 



r „„„2 



cos a 



(20) 



ahl + y/1 + (1 - cfi)hl 

When the surface is horizontal, h r = 0, we have / = 1, and II reduces to the expression 
(HE). 

A peculiar property of the expression (j2*U|) facilitating practical calculations is that it is 
quite insensitive with respect to variations in the polarization angle a, especially in the case 
when a is close to unity, which is in practice most important. Thus if we draw curves for 
H(9i) versus 9i for various input values of a in the whole region < a < 90° (not shown 
here) , we will find that the curves lie close to each other. For practical calculations involving 
unpolarized light it is thus legitimate to replace f(h r , a) by its average with respect to a. 
As (sin 2 a) = (cos 2 a) = 1/2, we can then write the surface force density as 

2n\I(r) 1 — a 



n 



-f(h r )n, 



c 1 + a 

where f(h r ) is equal to f(a, h r ) averaged over a, 

1 + (2 - a 2 )hl + ah 2 ^/l + (l-a 2 )h 2 + h 



(21) 



f(hr. 



1 + a 



.1 2 



a + + a 2 )h 2 ah 2 + y/l + (1 - a 2 )k. 



r\ 2 ' 



(22) 



This expression is valid also in the case of hydrodynamic nonlinearity. Note again that 
II is the normally-directed force per unit area of the oblique liquid surface. 

Finally, let us consider the force balance for the liquid column, assuming stationary 
conditions. When n 2 > n\ the surface tension force which acts upward, has to balance 
the combined effect of gravity and electromagnetic surface force, which both act downward. 
When the surface is given as h = h(r, 6), the mean curvature can be written as 



1 1 1 



rh r 



+ 



1 d 



ha 



Ri R2 rdr ^i + h 2 r + (h e /r) 2 r 2 36 + h 2 + {hg/r) 



(23) 



with sign conventions the same as in Eq. (|5]l. Thus for azimuthal symmetry, 

1 d rh r 
r 



1 1 

R\ R2 



and the force balance becomes 



26 



31] 



dr^TTh 2 ' 



(24) 



(Pi - P2)gh{r) - 



a d 
r dr 



rh T 



VT+h 2 



U(r). 



(25) 
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This equation follows from considering the equilibrium of a liquid column having unit base 
area. Here U(r) is the pressure scalar, i.e., n(r) = Il(r)n. Thus U(r) < 0. 

What expression to insert for U(r) in Eq. (J25|) . depends on the physical circumstances. 
Thus in the case of an unpolarized laser beam, we may use either the expression (|14j) with 
R = (R), T = (T), or alternatively use the expression (|21j). We will follow the latter option 
here. As noted, there is no restriction imposed on the magnitude of the slope of the surface. 



III. SOLUTION OF THE NONLINEAR EQUATION 

It is advantageous to introduce nondimensional variables. Let us first define the capillary 
length lc and the Bond number B, 

fc = ,/=^= B=OV, (26) 



(pi - pi)g \lc 

uj being the radius of the beam waist. The Bond number describes the strength of buoy- 
ancy relative to the Laplace force. If B <C 1, gravity is much weaker than the Laplace 
force. (The Bordeaux experiments covered the region 1CT 3 < B < 4.) We then define the 
nondimensional radius R and the nondimensional height H(R) as 

R=~, H(R) = ^1. (27) 

The fact that in practice a = ni/n 2 is very close to one, makes it at first sight possible 
to simplify the right hand side of the governing equation (}25|) . Namely, from Eq. (}22|) one 
would expect that f(h r ) — > 1. However, the situation is more delicate due to nonlinearity: 
if we keep f(h r ) in the formalism and calculate the elevation h(r), we will find that / gets 
a pronounced dip within the region where the beam is strong. Typically, if we draw a curve 
for / = f(R) versus the nondimensional radius R, we will see that / drops from 1 to about 
0.3 when R lies about 0.5. Mathematically, this is because the high steepness of the surface 
makes h r (or Hr) large enough to influence / significantly in a narrow region even when a 
is close to unity. Assuming a Gaussian incident beam, 

I(r) = ^e-**M, (28) 

7TUJ 

with P the beam power, we may write the governing equation (|23|) as a nonlinear differential 
equation for H: 

T t RR -f ^tl R -f ^rU r j-i„-2R 2 



Hrr + r-Hr + s^pH 3 
(1 + i**) 3 ' 2 



BH - — R — Mp* = —Fe f(Hfi). (29) 
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Here F is a positive constant at fixed temperature, 

2(n 2 -m)P 



^cg(px -P2)lc' 

which can for practical purposes be written as, since (dn/dp)T = —1.22 x 10~ 4 m 3 /kg 



(30) 



F _ 2 / dn\ P _ 26AQQ P (31) 
ircg V dp) T l 3 c l 3 c 

In the last equation, the dimension of lc is pm. 

The two quantities lc and a will vary with the temperature T in accordance with the 
theory of critical phenomena. Thus for the density contrast Ap = p\ — pi we have 



T-T c 



P 



Ap = (Ap) Q (— , (32) 

where (3 = 0.325, (Ap)o = 285 kg/m 3 , Tc =308.15 K being the critical temperature above 
which the mixture separates into two different phases. Similarly 



T-T c 



2;/ 



a = a ° \ T c ) ' ^ 
with v = 0.63, a = 1.04 x 10^ 4 N/m. More details can be found in Refs. HQ- We give 
here the practically useful formulas for a and lc- 

(T T \ ^ ,3 ^^ 
' (34) 

/rp _rp \ 0.468 

lc = 193 y-^jr^) M- ( 35 ) 

These two quantities are the only parameters that vary with temperature. There are thus 
three parameters in all in the problem, namely T, the beam power P, and the beam waist 
loq. Nondimensionally, the last two parameters correspond to F and B, Eqs. (13 1|) and (|26j). 

A. Numerical solution 

It is convenient to let H be positive downwards, so that in the formalism below we will 
replace H with —H. 

We start from the nondimensional governing equation in the form 

1 d I RH fi \ .>/>•• 



RdR 1^1 + H 2 R /B 



BH = —Fe~ f(Hji), (36) 
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with boundary conditions 

H R (0) = 0, H(oo) = 0. (37) 
We use a two-point method to solve the nonlinear differential equation iteratively. Define 



K = Jl + H%/B, S = Fe' 2 *f{H R ), (38) 



and let G = H R /K. We obtain the following first order system: 



with boundary conditions 



§ = KG, (39) 



G(0) = 0, H{oo) = 0. (41) 



We linearize the equations by means of lagging, i.e., we use values for H R from the last 
iteration in the nonlinear functions K and S. Using a nonuniform grid with n grid points, 
we integrate the equations between two grid points j and j + 1, letting m be the midpoint 
and ARj the distance between the points. We obtain 

AR 

H j+1 - H 3 = K m —^-(Gj + G j+1 ), (42) 

IAD A D 

G j+1 - Gj + — -^-(Gj + G j+1 ) - B^{H 3 + H j+1 ) = -ARj S m . (43) 
Here H R in K m and 5 m are evaluated as 

where the ff's are values from the previous iteration. With n grid points there are n — 1 
intervals and n — 1 sets of equations. This confirms with the fact that there are 2n quantities 
H and G; since there are two boundary conditions there remain 2n — 2 equations in all. 

To start the iterations we give initial values for H = const. e~ R2 . To deal with the 
boundary condition at infinity, we need in practice to replace "infinity" with a finite upper 
limit R = R e . The solution falls off quite slowly with R, so to use the naive condition 
H(R e ) = would require R e to be very large. To avoid calculating the long tail of the 
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solution, we can find a better boundary condition by using the lowest order term in an 
asymptotic expansion for H. When R is large, S ~ e~ 2R2 , and H\ is very small so that 
S^O, y/l + H 2 R /B « 1. Equation (jU| becomes 

^ {RHr) - BH = - ( 45 ) 
To lowest order this equation has the asymptotic solution H ~ R~ 1 l 2 e~^~ B ~ R ', which in turn 
implies that 

G = H R = - ( ^— + Vb\ H. (46) 



N 2i2 

We take this condition to replace the condition H = at R = R e . 

We solved the discretized equations using a Block-Bidiagonal-Matrix-Algorithm, devel- 
oped by one of the authors (S.H.). Our programming language was MATLAB. 

B. Results 

First, the following question naturally arises: at which powers P will the nonlinear cor- 
rection begin to be important? And what magnitudes of the centerline deformation (r = 0) 
does this correspond to? To get insight into this issue we constructed a number of figures 
(not shown here) for the surface height h(r) versus r, for various temperature differences 
T — Tc and beam radii ujq, and for various values of P. For each parameter set we made the 
calculation in two ways, viz. when taking the nonlinear correction into account, and when 
omitting it. Of course, there is a transitional region and no sharp limit distinguishing the 
linear and nonlinear regions. Our conclusions, based upon visual inspection of the curves, 
were that under normal conditions the linear region can be taken to prevail until P ~ 200 
mW. When P > 300 mW, nonlinear effects turn up. Generally, the nonlinear deformations 
are higher than the linear ones. (To give an example: at P = 220 mW the nonlinear cen- 
terline deformation was found to be only 2% higher than the linear deformation, whereas at 
P = 330 mW it was 5% higher .) 

we recall that in the linear regime we have the following simple formula 



From Refs. 
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for the centerline height h (r) (here in physical variables): 

ft ""« = 2^(S) T frT^^* 8 ^- (47) 

In the following we show three figures, each of them corresponding to given values of T — Tc 
and u . Each figure is based upon a full nonlinear calculation. First, Fig. 2 shows how 
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h(r) varies with r when T — Tc = 2.5 K and u>q = 4.8 /im. According to Eqs. (|35|). ()26|) 
this corresponds to Z c =20.3 /im, 5=0.0560. The powers are P = {300,600, 1200} mW. As 
mentioned above, we had to replace "infinity" with a finite outer limit R e . Numerical trials 
observing the sensitivity of calculated centerline deformations showed that the choice R e = 9 
was large enough. (For instance, an increase of R e from 9 to 10 would lead only to minute 
differences, the first three digits in the centerline deformation being the same.) Because of 
the cylindrical symmetry, only one half of the displacements (r > in the figure) need to 
be shown. It is seen that both powers 600 mW and 1200 mW lead to deformations much 
greater than 20 /im, and are clearly in the nonlinear region. The theoretical deflections 
for the three given values of the power are read off from the figure to be {12, 56, 112} /im, 
respectively. 

Our choice of input parameters makes the figure directly comparable to Fig. 6.1 in Cas- 
ner's thesis [24]. The experimental centerline displacements estimated from the photos are 
about 10 /im for P = 300 mW, 45 /im for P = 600 mW, and 75 /im for P = 1200 mW. The 
theoretical predictions are thus in this case larger than the observed ones, especially for the 
highest value of P. Moreover, one difference which is most noticeable is the absence of the 
observed " shoulder" in the theoretical solution in the case of large P. The shoulder occurs 
experimentally when the laser illumination is from below. There is at present no theory 
capable of describing this phenomenon. Mathematically, the establishment of the shoulder 
seems to be related to an instability; the real deflection jumps from one class of solutions 
of the nonlinear differential equation to another class. Video records actually show "jumps" 
in the surface when it gets formed, thus supporting our conjecture about an instability phe- 
nomenon. As for the observed width of the surface displacement, there is good agreement 
with the theoretical prediction. 

Fig. 3 shows analogous results for the case T — T c = 2.5 K, u = 8.9 /im, thus a consid- 
erably broader beam. Here lc =20.3 /im, B = 0.193. The powers are P = {360,600,890} 
mW. In this case, the value R e = 7 was found to be sufficient, for the same reasons as above. 
Our results c m b e _ d with Fl , ,3 <W thesis Q. The tWtica! centerhne 
displacements for the three mentioned cases of P are {10, 19,47} /im, which all agree well 
with the observed values. Also in this case there occurs a shoulder experimentally, but it is 
not so pronounced as in the previous case. 

Finally, in Fig. 4 we show the case T — Tc = 3 K, ujq = 5.3 /im, corresponding to lc— 22.1 
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fim, B = 0.0576, for powers P = {300, 590, 830} mW. Again, the outer nondimensional 
radius R e = 9 was found to be appropriate. The theoretical centerline deflections are now 
seen from the figure to be {10,39,65} /im, respectively, for the given values of P. We 

l: the 



may compare this with the photos shown in Fig. 2 in Ref. [27] or Fig. VI. 5 in 
corresponding experimental dispacements are about {10,40,55} /im. Also this time we see 
that the agreement between theory and experiment is good for low powers, but that the 
theory over-predicts the displacement when the power gets large. The last-mentioned effect 
is generally most pronounced when the radius of the laser beam is small. 



IV. SUMMARY, AND FINAL REMARKS 

The "giant" deformations of fluid interfaces recently obtained in the experiments of Cas- 



ner and Delville 
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24 



25 



27. 



281 ] with the use of moderate laser beam powers P (P 
typically lying between 500 and 1000 mW) have accentuated the need of taking into account 
nonlinear geometrical effects in the theoretical description of the interface deformation. As a 
rule of thumb, inferred from a visual inspection of the figures, nonlinear effects are expected 
to be appreciable when the deformations become larger than about 15 /im. When the radius 
uj of the laser beam is small, typically u ~ 5 //m, a power P ~ 1 watt may easily lead to 
deflections around 100 /xm. The physical reason for the occurrence of giant deformations 
is the lowering of surface tension caused by working with a phase-separated liquid mixture 
close to the critical point. 

The nonlinear governing equation in nondimensional form can be taken as in Eq. (|2*9l 
or, what was found more convenient, as in Eq. (|36|) where G = Hr/K is considered as the 
dependent variable. We solved the set of equations (j3~9"j) - (j4*Tj) numerically. As a consistency 
check, we made also an analogous calculation starting from Eq. (|29|) . and got the same 
results. Figures 2-4 show some examples of our calculations; these are all directly comparable 
with the Casner-Delville experiment. 

Some general conclusions that can be made from our calculations, are the following: 

1. For given values of T — Tq and ujq, a larger P causes a larger deformation. 

2. For a given T — Tc, a smaller ujq causes a larger and narrower deformation. 

3. For a given ujq, a smaller T — Tc causes a larger and narrower deformation. 

4. Very large beam waists (a>o ~ 20—30 /im) are not able to cause a nonlinear deformation, 
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not even for the largest P and smallest T—Tc investigated in the Casner-Delville experiment. 

5. For small T — Tc and small uq (for instance T — Tc = 2.5 K and coq = 4.5 /mi), a power 
of 300 mW is not enough to cause a nonlinear deformation. However, a further decrease 
in temperature, such as to the value T — Tc=1.5 K, will take also the 300 mW-induced 
deformation into the nonlinear regime. 

All the items listed above are expected on physical grounds. A large incident power 
concentrated on a narrow cross section means a large electromagnetic field intensity and 
thus a large surface force. The enhanced deformation for small T — Tq is due to the fact 
that the restoring buoyancy force (~ Ap) and Laplace force (~ a) vanish at T = Tc ; cf. 
Eqs. flsgD and (jHHj). 

Concretely, when comparing our results with the Casner-Delville observations, we find 
that for broad beams the agreement between theory and experiment is quite good; cf. our 
discussion of Fig. 3 above. There is however a considerable theoretical over-prediction of the 
deflection in the case of narrow beams and high powers. Most strikingly, this is shown in the 
case a>o=4.8 /im, P=1200 mW, as discussed in connection with Fig. 2. The physical reason 
for this discrepancy is not known. It may be related to the production of heat in the liquid 
in the presence of the strong field, or to the loss of radiation energy because of scattering 
from the non-avoidable corrugations on the liquid interface. Perhaps the most intriguing 
possibility is that the discrepancy is related to the reflection of radiation energy from the 
interface "shoulder", which is seen to be produced in strong fields when the illumination is 
from below. This effect is most likely related to an instability; the system decides to switch 
from one class of solutions of the nonlinear governing equation to another class. To our 
knowledge, no explanation exists of this effect. 
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Figures 
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FIG. 1: Definition sketch of the displaced surface. The laser illumination is from below. 




r(um) 

FIG. 2: Theoretical height h(r) of displaced surface versus radius r when T — Tq = 2.5 K, ujq = 
yum, for three different laser powers P. The undisturbed surface is at h(r) = 0. 
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FIG. 3: Same as Fig. 2, but with T - T c = 2.5 K, uj = 8.9 fim. 
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FIG. 4: Same as Fig. 2, but with T -T c = 3K,lo = 5.3 /im. 
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